mod returns the floating point remainder of the
division of the first argument by the second argument.
If the divisor is zero or if the division would overflow, then it returns
zero. When the arguments mod are both matrices, then an
element by element operation is performed, so they must have the
same dimensions.